Density Matrix From Photon Number Tomography 



Stefano Mancini and Paolo Tombesi 

Dipartimento di Matematica e Fisica, 
Universitd di Camerino, 1-62032 Camerino, Italy 

and 

Istituto Nazionale di Fisica della Materia, Italy 



Vladimir 1. Man'ko [] 
Osservatorio Astronomico di Capodimonte, 1-80131 Napoli, Italy 



Abstract 

We provide a simple analytic relation which connects the density operator 
of the radiation field with the number probabilities. The problem of experimen- 
tally "sampling" a general matrix elements is studied, and the deleterious effects 
of nonunit quantum efficiency in the detection process are analyzed showing 
how they can be reduced by using the squeezing technique. The obtained result 
is particulary useful for intracavity field reconstruction states. 
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After the seminal paper of Vogel and Risken [|1| and the experimental realization 
of their result by Smithey et al. |^ it becames clear that the homodyne measurement 
of an electromagnetic field permits the reconstruction of the Wigner function of a 
quantum state, by just varying the phase of the local oscillator, since then named 
optical homodyne tomography. In Ref. it was shown, indeed, that the rotated 
quadrature distribution may be expressed in terms of any s-parametrized Wigner 
function 0. 

These papers give rise to a pletora of other papers |^ where the limitations due 
to the detectors quantum efficiency r] were also taken into account, and a fundamental 
limit was established showing the impossibility of the tomographic reconstruction for 
detectors quantum efficiency less than 0.5 0. Very recently it was shown |^ that, 
in principle, there is no such limit; however, in order to get the desired results for r] 
arbitrarily small, one needs to employ a rather complicate mathematical procedure 
for loss error compensation, which in practice make still persistent a lower limit on r] 
values because of numerical problems. Anyway, all the above mentioned papers rely 
on homodyne measurements, or any other phase dependent measurement process. 
While this work was in progress we became aware that Wallentowitz and Vogel [§] 
proposed the s-parametrized Wigner function reconstruction similar, in its essence, 
to the present one by using direct photon counting and contemporarely an analogous 
scheme was adopted by Banaszek and Wodkiewicz |P; however, in the following we 
shall discuss the relevant differences with our work. 

We shall show in this letter that the state retrieval by photon counting can be used 
not only when the output beam is mixed with a reference field at a beamsplitter, as 
m Refs. but also in a physical situation similar to the one proposed by Brune et 

al. [|lOl for the generation and measurement of a Schrodinger cat state allowing thus 
the possibility of reconstructing its density matrix and, more generally, the possibility 
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to reconstruct cavity QED field state when there is not a direct access to the field. 

We shall show that the sampling of density matrix elements can be achieved by 
photon number measurements, for detectors quantum efficiency greater than 0.5, 
in any base but coordinate representation. This rj > 0.5 limitation, however, can be 
beaten if we use a physical technique instead of applying the mathematical procedure 
ofRef. §. 

We shall call the procedure Photon Number Tomography because it permits the 
reconstruction of the density matrix elements, in a suitable representation, by just 
detecting the number of photons at the given reference field and then scanning both 
its phase and its amplitude; differently from the usual homodyne tomography where 
a marginal distribution is recorded by homodyne measurements and then scanning 
only the phase. 

As claimed in Ref. 0, it is possible to write 

(Pa, 



P 



f—W{a,s)f{a,-s), (1) 
J vr 

where the s-ordered wheight function W{a,s) may be identified with the quasiprob- 
ability distributions Q{a), W{a) and P{a) when the ordering parameter s assumes 
the values —1, 0, 1 respectively; while the operator T represents the complex Fourier 
transform of the s-ordered displacement operator D{^,s) = l)(^)e*l^l^/^, which can 



also be written as 11 



T{a, s) = Y^^Dia) j D-\a) . (2) 

On the other hand the weight function s) is the expectation value of the operator 
T{a, s) 0, i.e. W{a, s) = Tr{pT(a, s)}, then we obtain 

P = / — E 1 [-^) {n\D{a)pD~\a)\n)T{-a, -s) . (3) 

J vr 1 - s Vs - 1/ 

Thus, the weight function W{—a, s) is related to the ability of measuring the quantity 
{n\D{a)pD~^{a)\n) by scanning the whole phase space |]12|, i.e. by just varying 



a. Now, we may consider one mode of the radiation, whose state p one wants to 
reconstruct, contained inside a cavity and, immediately before the photon number 
measurement, a coherent reference field is "added" |13|, so that we may recognize 

P(n,a) = Tr{Z)(a)pZ)-^(«)|n)(n|} = (n|Z)(a)pD-^(a)|n) (4) 

as the probability to detect n photons after the injection of the reference field a. 
The addition process we are considering, following Ref. |TD[, "is quite different from 



the combination of fields produced by a beam splitter, which mixes together distinct 
modes coupled to its two ports and introduces vacuum noise even in the absence 
of any classical input field". We are indeed describing a much simpler field ampli- 
tude superposition mechanism, discussed in the Glauber's pioneering work |13| . The 
photon number distribution (|) results as the projection of the field state p over a 
displaced number state . 



The photon counting could be made either by means of atoms [jT5| as in the case 
of microwave cavity field or by direct detection of the outgoing optical total field. 
Furthermore, setting 

K,{n,(x) = - T(-a,-s), (5) 

\ — s Vs — 1/ 

we have, from Eq. (H) 

°° r (P'cy 

p=Y. —nn,a)Ksin,a). (6) 

Thus, analogously to Ref. [Q, we may assert that a density matrix element can be 
experimentally sampled if there exist at least one value of the parameter s inside the 
range [—1,1] for which the corresponding matrix element of the kernel operator Kg 
is bounded. From Eqs. (||) and (H), one immediately recognizes that this is possible 
for s G (—1,0], in the number, coherent and squeezed representations and not in 
the coordinate basis. For s = —1, in Eq. (^) n = only survives, however as was 



shown in Ref. 0, the operator T becomes singular and can be used to construct an 
arbitrary density matrix when is only weighted with the well behavied function Q. 
It means that V{0, a) = Q{a) in Eq. (P). 

Let us now consider the more realistic case of nonunit quantum efficiency rj. Ac- 



cordingly to Ref. ||T^ we have 

k 

V 



where represents the same of (^) but in presence of efficiency rj, i.e. 

P^(n,a) = Tr jD(a)p D-\a) : 6"""''^^^^^^ :| . (8) 



ni 



If we substitute Eq. into Eq. we obtain 



m=0 



TT 



1=0 



r] 



(9) 



where the quantity inside the square brackets can be considered as a modified kernel 

K,^^{m,a) = -^( '^''~'^^^ ] f{-a,-s). (10) 
1 — s \ rjs — 7] J 

Again, this kernel results bounded only in the number, coherent and squeezed rep- 
resentations, but in that case the s values for which that occours are determined by 
the quantum efficiency; in particular we have from Eq. (p!0| ) s G (— 1, — (1 — t])///], so 
that 7] should be greater than 0.5, similar to Ref. 0. 

It is also interesting to note as was stressed in Ref. p, that if one chooses 
s = l — 2/r7in Eq. (pUj), only the term with m = survives in Eq. (y), and this 
means that the desired density matrix elements are simply related to the zero-count 
probability which gives only the Q-function. Unfortunately, only quantum efficiency 
f] = 1 (i.e. s = —1) allows to get the Q-function in this simpler manner, otherwise, 
for 1] < 1, s takes forbidden values (i.e. becomes smaller than —1). 
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Let us now show how the above hmit on detectors efficiency can be beaten. In 
case one squeezes the radiation inside the cavity immediately before the injection of 
the reference field, the superposition mechanism O] can be read as 



D{a)SiC)pS-\C)D-\a) = b{a)pb-\a) , (11) 



where S{() is the well known squeeze operator with the squeezing parameter 
C = ICk"''- Of course one has to be sure that the squeezing process is like a kick during 
which the natural evolution of the system is negligible; only under this assumption 
Eq. ( [TlD holds. It could be physically realized, for example, by using a (5-kicking of 
frequency of the cavity mode [|T^ . 

Since the weight function W{a,s) can be written as the complex Fourier transform 
of the characteristic function, let us examine the relation between the operator of Eq. 
([Tl| ) and b{a)f)b~^{a) in terms of characteristic functions. By definition we have 

X(e, s) = Tr{p~D(e, s)} = i:i{pS-\C)e^''e-^''S{C)}e-\^\""e^\^\'" . (12) 

Remembering how the squeeze operator acts on the annihilation (creation) operator 



T7| , we obtain 

X(e, s) = Tiipmfi* + eV)}e^l«l'/' = + eV)e^l^l'/' , (13) 

with fi = cosh \ (\ and u = e*''' sinh \ (\. Writing now the variable ^ in polar coordinate 
^ = l^le*"^ and locking the squeezing parameter phase with that of the reference field, 
i.e. ip = 20, we get ^/i* + ^*z/ = ^(cosh|C| + sinh|C|) = ^e'^' = ^A, such that 
X{^, s) = x(^/A, sA2) and W{a, s) = A'^W{Aa, sA^). Hence, from Eqs. (|lD, ©, (0) 
and the last expression, we get 

°° r (Pa 

P = E / ^vi^^ a)ICs,r,in, a/A) , (14) 

where 

A^ = 1 ,sA^-, ) A ' ' ^''^ 



from which one immediately recognizes that the ordering parameter s in the power 
term, is now scaled by the factor with respect to the one in Eq. ([To|) , thus it 
ranges in the interval s G (—1, —(1 — r])/riA'^], allowing values of quantum efficiency 
lower than 0.5. It is then possible to define an effective quantum efficiency fi = 
A^/[A^ + (1 — r])/r]] which could be close to unit for sufficientely large squeezing 
independently of the real quantum efficiency. 

For an easier implementation of this scheme, we could choose in Eq. ([I5D s = 
A^^(l — 2/1]), in order to eliminate the sum over n in Eq. (0), directly relating the 
density matrix elements in any allowed basis to the zero-count probability; but in this 
case, due to the presence of the factor A~^, the allowed values of s are not s = — 1 as 
in Ref. ||^, which only permits the reconstruction of the more smoothed Q-function, 
rather the value of s can approach zero from below depending on the value of the 
squeezing parameter, independently of rj. 

Thus, by squeezing preventively the field in a cavity allows one to detect the state 
of radiation inside it even for low efficiency detection. The use of squeezing technique 
to reduce the effect of nonunit quantum efficiency was already proposed in Ref. f^, 
but we use it in a different way. In Ref. [^, indeed, antisqueezing preamplification 
of the signal at the beam splitter was proposed rather than the antisqueezed signal 
inside the cavity. We would stress the fact that our Photon Number Tomography 
scheme becomes especially useful in the intracavity optical tomography where other 
similar schemes p, |^ are not applicable; in fact it can be adopted in a situation in 
which the photon number is measured indirectly using a "sequence" of atoms passing 



through the cavity [O] , with the quantum efficiency, determined only by the duration 
of the measurement process (i.e. the length of the "sequence"), that could be very 
high [|l^]. In that case, the scheme has also the advantage of being QND. Finally our 
scheme results suitable for cavity QED characterization, like Ref. [|T9l, allowing the 
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reconstruction of nonclassical states as well, which are extremely sensitive to losses 
and then their detection seem prohibitive by means of an outgoing field as in Refs. 

At last, our approach could be considered an extension of the known mathemati- 
cal principle of tomography. In fact given a density operator p and a group element 

one can create different types of tomography if, by knowing the matrix elements 
{x\Qf)G~^\x) from measurements, is able to invert the formula expressing the den- 
sity operator in terms of the above distribution (the x may denote either continous 
or discrete eigenvalues). For the inversion procedure one can use the properties of 
summation or integration over group parameters Q. The known tomographies are 
just given by this construction; for Radon transform or homodyne tomography, Q is 
the rotation group, for symplectic tomography |^ it is the symplectic group, but in 



principle one can use other groups. The only problem is mathematical one to make 
the inversion and/or physical one to realize the transformation Q in laboratory. 
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